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\Q ' We study a configuration of D-branes in string theory that is described at low energies by 

a four-dimensional field theory with a dynamically broken chiral symmetry. In a certain 
region of the parameter space of the brane configuration the low-energy theory is a non- 
g_) ' local generalization of the Nambu-Jona-Lasinio (NJL) model. This vector model is exactly 
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X, 

H , D-branes living in the near-horizon geometry of other branes. In a different region of 

parameter space the brane construction gives rise to large Nc QCD. Thus the D-brane 
system interpolates between NJL and QCD. 



solvable at large A''^ and dynamically breaks chiral symmetry at arbitrarily weak 't Hooft 
coupling. At strong coupling the dynamics is determined by the low-energy theory on 
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1. Introduction 

Quantum chromodynamics (QCD), the theory of the strong interactions, is weakly 
coupled at high energies but strongly coupled at the scale of typical hadron masses (~ 1 
GeV). It has proven difficult to use this theory to study analytically the properties of 
low- lying mesons and baryons. 

In the approximation in which Nf ffavors of quarks are taken to be massless, the 
Lagrangian of QCD has a chiral global symmetry U{Nf)L x U{Nf)ji, acting on the left 
and right-handed quarks, ql, Qr- This symmetry is expected to be dynamically broken to 
the diagonal subgroup t/(A^/)diag- Analyzing this breaking is difficult due to the strongly 
coupled nature of the theory. 

An important early example of dynamical chiral symmetry breaking occurs in the 
Nambu-Jona-Lasinio (NJL) model [|l|. This model contains fermions ql, qr, which trans- 
form in the fundamental representation of a global U{Nc) symmetry, and interact via a 
local four Fermi coupling, 

>Cint = Gq[ ■ qnq]^ ■ Ql , (1.1) 

where the dot product denotes a contraction of the color indices. The fermions ql and qr 
also transform in the fundamental representation of U^Nf)^ and U{Nf)fi symmetries of 
£int, respectively. 

The four Fermi interaction (|1 . 1|) provides an attractive force between the left and 
right-handed quarks. The idea of |jl[], based on an analogy with the BCS treatment of 
superconductivity, was that this force may destabilize the trivial vacuum and lead to the 
breaking of chiral symmetry. To test whether this indeed occurs in the NJL model one 



must work with a finite UV cut-off, since the interaction ( |1 . 1| ) is non-renormalizable. In 
^ it was shown that the resulting model breaks chiral symmetry when the coupling G 
exceeds a certain critical value. 

The non-renormalizability of the NJL model implies that many of its predictions are 
sensitive to the precise nature of the UV cut-off. Thus, it is natural to look for a renor- 
malizable ffeld theory in which the ideas of [|l| can be tested in a more controlled setting. 
Such a setting is the two-dimensional analog of the NJL model, which is asymptotically 
free. Moreover, since it is a vector modelu it is exactly solvable at large N,. and ffnite Nf. 
This model was analyzed by D. Gross and A. Neveu [Q, who showed that it breaks chiral 
symmetry and generates a mass scale via dimensional transmutation. 



See e.g. H for a review of large A^ vector models. 



The Gross- Neveu model provides a beautiful confirmation of the ideas of Nambu and 
Jona-Lasinio, and makes it interesting to look for more realistic four-dimensional models 
with similar properties. Restricting to asymptotically free theories one is naturally led to 
QCD, since there are no asymptotically free field theories in four dimensions without gauge 
fields. Unfortunately QCD, like most models that involve Nc x Nc matrices, is difficult to 
solve even in the limit Nc —>■ oo. 

In this paper we explore a different route. We study a configuration of D-branes in 
string theory that reduces at low energies to a field theory of left and right-handed fermions 
Ql, Qr, with an adjustable attractive interaction. The brane configuration we start with 
preserves a U{Nf)L x U{Nf)^ global symmetry but, as we will see, this configuration 
(or vacuum) is unstable for all values of the coupling. The true vacuum has a non-zero 
condensate {q\ ■ qr) . 

At weak coupling, the infrared dynamics of ql, Qr is captured by a non-local NJL 
model - a vector model that does not contain dynamical gauge fields and can be solved 
exactly in the limit Nc —>■ oo, like the Gross-Neveu model. It breaks chiral symmetry 
for arbitrarily weak coupling and generates a mass for ql, qr. At strong coupling, the 
useful description is in terms of D-brane dynamics in curved spacetime, which can also be 
analyzed quite explicitly and exhibits similar properties. 

The embedding of the NJL model in string theory that we describe provides another 
bonus. Despite its non-renormalizability, the NJL model (|1.1|) and its non-local general- 



izations have been used extensively in phenomenological studies of hadrons, and appear to 
give a rather accurate description of their properties (see e.g. [^|] for reviews and further 
references). However, it is not understood to what degree these models can be thought of 
as effective low-energy descriptions of QCD. The realization of NJL as a low-energy limit 
of the theory on D-branes provides a new perspective on this problem. By varying the 
parameters of the D-brane configuration, one can interpolate between NJL and QCD. It 
is likely that the interpolation is smooth, i.e. the two models are in the same universality 
class. 

The idea that embedding field theories in string theory makes it easier to under- 
stand otherwise mysterious phenomena in field theory is familiar from other contexts. For 
example, thinking of A^ = 4 SYM in four dimensions as the low-energy limit of the six di- 
mensional (2, 0) SCFT compactified on a two-torus makes manifest the SL{2, Z) S-duality 
symmetry of A = 4 SYM, which acts as the modular group of the torus. Thinking of 



A^ = 1,2 SYM as low-energy limits of the (2,0) theory wrapped on other Riemann sur- 
faces provides a geometric realization of the Seiberg duality of A^ = 1 SYM and of the 
Seiberg-Witten curves of A^ = 2 SYM (see e.g. for a review). Similarly, in our case 
both NJL and QCD are realized as low-energy limits of the compactified (2, 0) theory in 
the presence of defects. It is interesting that the (2, 0) theory plays a central role in all of 
these constructions. 

The plan of the paper is the following. In section 2 we describe the brane configuration 
and its massless excitations. We introduce the fundamental length scales in the problem 
and discuss the range of validity of the weak and strong-coupling approximations. In 
section 3 we study the low-energy dynamics of the branes at weak coupling and show 
that it is described by a non-local NJL model which exhibits chiral symmetry breaking at 
arbitrarily weak coupling (for large A''^). 

In section 4 we study the low-energy dynamics at strong coupling using a description 
in terms of D-branes in curved spacetime. We find the chiral condensate and comment on 
the relation to the weak coupling analysis. In section 5 we discuss a generalization of the 
D-brane construction which interpolates between NJL and QCD. We conclude in section 6 
with comments on our results and a discussion of possible extensions of our analysis. Our 
conventions and a few results used in the text are described in an appendix. 

2. The D-brane configuration and some of its properties 

Motivated by earlier studies of string theory duals of QCD in the large Nc limit with 
the number of flavors Nf held fixed (see e.g. [1T0| - |T3[] ) , we will consider a brane configuration 
in IR^'^ that includes three kinds of D-branes: DA, D8 and D8. The different branes are 
extended in the directions 

0123456789 

D4 : X X X X X (2.1) 

-D8, D8 :xxxx xxxxx 

and are arrangedcl as indicated in figure 1. Classically, the Nf D8 and D8-branes are 
parallel and separated by a distance L in the x^ direction. They intersect the Nc -D4- 
branes at the origin of the IR labeled by (x^, ■ ■ ■ , x^) and at x^ = ±L/2. We will see that 
this classical picture is modified by quantum effects. 



This brane configuration was studied in |13], with the x direction compact and with anti- 
periodic boundary conditions for the ferniions. As we will see, the hmit where the radius of x^ 
goes to infinity is particularly instructive. We will return to the compact case in section 5. 
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Fig. 1: The configuration of Nc D4-branes and Nj D8 and DS-hvanes wliich leads 
to QCD in one limit and the NJL model in another. 

Since the different branes share the directions (x*', x^, x^, x"^), it is natural to study 
the low-energy dynamics seen by a 3 + 1-dimensional observer living in the intersection 
region. In the following sections we will do this in different regimes of the parameter space 
of the brane configuration. Here, as preparation for this study, we start with some general 
observations on this problem. 

As usual in D-brane physics, there are different sectors of open strings that we need 
to consider, p — p strings {p = 4, 8, 8) give rise to p + 1-dimensional gauge theories with 
gauge groups U{Nc), U{Nf)L and U{Nf)ii respectively. The field content of these gauge 
theories is obtained by dimensional reduction of A^ = 1 SYM in ten dimensions. The gauge 
couplings are (see e.g. [jl^) 

{2nr-^gjr' ■ (2.2) 



2 

9p+i 



Since these theories live in more than 3-1-1 dimensions they are not of direct interest to 
us, as we are looking for modes localized in the extra dimensions. The 



f/(A,) X U{Nf)L X UiNf) 



R 



(2.3) 



gauge symmetry of the D-branes gives rise to a global symmetry of the 3 -|- 1-dimensional 
theory. 

Normalizable modes in 3-1-1 dimensions are found in the 4 — 8 and 4 — 8 sectors, which 
contain states localized near the corresponding brane intersection. Open strings in these 
sectors have six Dirichlet-Neumann directions. The NS sector states are massive; the only 
massless modes are spacetime fermions. 4 — 8 strings give a left-handed Weyl fermion q^ 



localized at x^ ~ —L/2, which transforms in the {N^^ A^/, 1) of the global symmetry ( |2.3|) , 
while 4 — 8 strings give a right-handed fermion localized at x^ ~ L/2, which transforms 
as (A^c, 1, Nf). A nice feature of the brane setup of |T^] is that the left and right-handed 
fermions are separated in the extra dimensions. 

The low-energy dynamics of the D4 — D8 — D8 system can be formulated in terms of 
ql and qji. Alternatively, one can integrate out these fields and study the dynamics of the 
gauge fields in the presence of a localized source (see [^ for a related recent discussion). 
We will follow the former approach, which is closer in spirit to the standard field theoretic 
treatment of 3 + 1-dimensional dynamics. 

Consider first the case of a single intersection, of Nc il>4-branes and Nf I?8-branes, 
which can be thought of as corresponding to the limit L -^ oo of the configuration of figure 
1. We would like to study this system in the limit Qs —^ 0, Nc — > c>o with QgNc and Nf 
held fixed, and at energies much below the string scale. In this limit the only interactions 
between the fermions ql are due to their non-derivative couplings to the Z)4-brane gauge 
fields. The leading effect of these couplings is a four-Fermi interaction whose strength is 
proportional to g'^, and whose precise form will be described in the next section. 

While the fermions ql live in 3 -|- 1 dimensions, the gauge fields that they exchange 
are 4 -|- 1-dimensional. The 't Hooft coupling of the five- dimensional gauge theory 

A = |iiV. , (2.4) 

has units of length. This length determines the range of the four-Fermi interaction. For 
example, one can show that this interaction gives rise to corrections to the propagator 
{q\{x)qL{y))-, which are suppressed relative to the leading, free field, result by powers of 
\/\x — y\. Thus, at distance scales much larger than A these interactions can be neglected, 
while in the opposite regime |a; — y| ^ A they are large. We will discuss both limits below. 

We next turn to the case where both the D8 and il>8-branes are present. Since we 
would like to study the low-energy theory on the branes, we will take the ratio of the 
distance between the eight-branes and the string length, L/ls, to be large but finite in the 
limit (7g — > 0. This will allow us to neglect effects such as the non-trivial dilaton and RR 
ten-form field strength created by the L)8-branes. 

For finite L there are interactions between the left-handed and right-handed quarks 
which are due to exchange of il>4-brane gauge bosons. These interactions are weak for small 
A/L, and become stronger as one increases this parameter. As in other brane systems. 



it is useful to label the different regimes by the value of the dimensionless parameter 
QsNc = X/ls- For small QgNc the hierarchy of scales in the problem is the following 

A < /, < L . (2.5) 

As mentioned above, we will only be interested in the physics at distance scales much larger 
than Is (but which may be larger or smaller than L) . In this regime we can neglect stringy 
effects on the dynamics of the fermions, as well as most of the effects of the dynamics of 
the gauge fields on the D4 and Z)8-branes. The only effect we need to keep is the one 
gluon exchange between the left and right-handed quarks, since it gives a force between 
them which is absent at infinite L. Unlike the propagator corrections mentioned above, 
this is not a small correction to an existing effect, but rather the leading interaction in this 
channel.0 

Thus, in the limit (|2.5| ) the dynamics is described by a field theory of the fermions ql 
and qn, with a non-local interaction due to one (five-dimensional) gluon exchange. The 
theory contains a natural UV cut-off Ig- Since this is much larger than A, the scale at 
which the non-linear dynamics of the five-dimensional gauge field becomes important, we 
do not have to include corrections to the Lagrangian obtained in the one gluon exchange 
approximation. In the next section we will study the resulting non-local NJL model, and 
show that its long distance dynamics is non-trivial. 

As QsNc increases, we get to the regime /^ ^ A ^ L, in which we can still treat 
ql and qn as weakly interacting via single gluon exchange, but the natural UV cut-off is 
now A and not Ig- Further increasing QsNc we get to A ~ L, where multi-gluon exchange 
processes can no longer be neglected, and the single gluon exchange approximation breaks 
down. For 

/, < L < A , (2.6) 

five-dimensional gauge theory effects give rise to a strongly attractive interaction between 
ql, Qr- This interaction can be studied using the results of [0, by analyzing the dynamics 
of DS-branes in the near-horizon geometry of the il>4-branes. This will be discussed in 
section 4. 



One can ask whether we should include additional interactions between ql and qr due to 
exchange of the scalars on the D4-branes. These interactions provide small corrections to those 
due to gluon exchange, since the scalars are derivatively coupled to the quarks. For the scalars 
$5, ■ • • , $9 this is due to the exact symmetry of the brane system $i — ;• $i + const. $4 is a 
component of the 4 + 1 dimensional gauge field, and is derivatively coupled to the quarks due to 
4 + 1 dimensional gauge invariance. 

6 



3. Weak coupling analysis 

In this section we study ttie D4 — D8 — D8 system in the weakly coupled regime ( |2.5|) . 
We start by deriving the low-energy effective action for the fermions ql, qn. Consider first 
the case of a single intersection, say that of A''^, il>4-branes and Nf DS-branes. Taking the 
intersection to be at x* = 0, the effective action for the Weyl fermion q^ and U{Nc) gauge 
field Am is given by 



S= (Fx 



45, 



l^Fl,^ + 5{x^)qla^{id^ + A^)qL 



(3.1) 



Note that while the first term (the gauge field Lagrangian) is integrated over the 4 + 1- 
dimensional worldvolume of the Z)4-branes with coordinates x^, M = 0,1,2,3,4, the 
fermion Lagrangian is restricted to the 3 -|- 1-dimensional intersection at x"^ = with 
coordinates x^, /U = 0, 1, 2, 3. 

Since the physical degree of freedom at the intersection is the fermion g^, it is useful 
to integrate out the five-dimensional gauge field Am- Keeping only the quadratic terms in 
the gauge field Lagrangian, and working in Feynman gauge, we find the following effective 
action for g^: 

Seff = i d'^xq^a^d^qi - ^^ / d'^xd!^yG{x - y, 0) qlixp^qLiy) q\{y)a^qL{x) 

(3.2) 
where G{x^, x'^) is proportional to the scalar propagator in 4 -f 1 dimensions, 

G{x'',x^) = i ^, (3.3) 

((a;^)^ — x^a;'^)2 



The color indices in (3^) are contracted in each term in brackets separately, while the 
flavor ones are contracted between q^ from one term and q]^ from the other. Thus the 
quantities in square brackets are color singlets and transform in the adjoint representation 
of U{Nf)L. To derive ( |3.2D we used a Fierz identity described in appendix A. 

The action ( [3?^ ) can be treated using standard large A^ techniques. The solution has 
the following structure. For distance scales much larger than A ( |2.4|) , the fleld q^ is free. 
The effects of the interaction grow as the distance scale decreases, and become important 
at distances of order A. 

At these distances there are two other types of interactions that are not taken into 
account in ( p.2| ). One is due to the non- linear terms in the gauge fleld Lagrangian ( |3.1| ). 



Their contributions are not small at distances of order A and adding them brings back the 
full complexity of the large Nc gauge theory. The other is due to string corrections. If the 
't Hooft parameter QgNc is small, the scale A ( |2.4| ) is much smaller than the string scale 
(see (^.51)). Thus, before we get to the point where the gauge theory effects become large, 
we reach a regime where the dynamics of ql is dominated by exchange of massive string 
states. In this section we will restrict the discussion to distances much larger than Ig and 
A, for which we can neglect all these effects, such that ql is free in this limit. 

We are now ready to discuss the case of interest, which contains D8 and L>8-branes 
separated by a distance L. The analog of the Lagrangian ( |3.1|) takes in this case the form 



S 



(Fx 



^niN + ^{^' + ^)cii-^^{^d^ + A^)qL + 5{x' - ^)q^Ra^{id, + A 



(3.4) 

Integrating out the gauge field in the single gluon exchange approximation and neglecting 
interactions of the form in (|3.2| ) we get 



Ses =i / d'^x [q[or'^d^qL + q^or'^d^qR 



a I r 1 r 

+ ^ / d'^xd'^yG{x - y, L) q\{x) ■ qniy) ql^{y) ■ qL^x) 



(3.5) 



where we again used a Fierz identity from appendix A. As mentioned above, ( p.5| ) provides 
an accurate description of the dynamics of q^-, qn when L and all other distance scales in 
the problem are much greater than A, /g. 

To solve the non-local Nambu-Jona-Lasinio model {^^ at large Nc it is convenient to 
introduce a complex scalar field T{x,y) which transforms as (1, Nf, Nf) under the global 
symmetry (^.3]), and rewrite the quartic interaction (|3.5|) as follows: 



Seff =i / d'^x (qlcr^d^qL + q]i<y^d^qR 
+ j d^xd^y 



A G{x — y, L) 

The equation of motion of T is 

A 



NcT{x,y)T{y,x) - | / ^ , t^/ \ U \ r \ 

+ T{y, x)ql{x) ■ qR{y) + T(x, y)qR{y) ■ qL{x) 



(3.6) 



(3.7) 



Plugging it back into ( p^ ) (or, equivalently, integrating out T, T), one recovers the origi- 
nal action (|3.5|). As usual in large A^ theories, we would like instead to integrate out the 



T{x, y) = j^G{x - y, L)q\{x) ■ qniy) 



fermions and obtain an effective action for the scalars, whose dynamics becomes classical 
in the limit Nc —>■ oo. Since we are mainly interested in the question of chiral symme- 
try breaking, we would like to compute the expectation value of T(x, y) in the vacuum. 
Poincare symmetry implies that the latter must be a function of (x — y)^. Thus, we will 
make this simplifying assumption and compute the effective action for T(x, y) = T{\x — y\). 
To leading order in the 1/Nc expansion, the only corrections to the classical action 
of T, T come from one loop diagrams with an arbitrary number of external T, T fields. 
Computing the one-loop effective potential in the standard way, adding to it the classical 
term from ( p.6|) and dropping an overall factor of Nc and the volume of spacetime leads 
to the effective potential 



. _^ _^ Jx2-FL2)3/2 f d^k , / T(k)T(k)\ 

where T{k) is the Fourier transform of T(x) (see appendix A for conventions). In ( |3.8| ) we 
have also Wick rotated to Euclidean spacetime. 

The equation of motion for T that follows from (|3.8| ) (the gap equation) is 

...(£!±i!)lr(.)e-^ . _|i|_ . (3.) 

The trivial solution of this equation, T{x) = 0, corresponds to a vacuum with vanishing 
{Ql ' ^r) ^^^ unbroken chiral symmetry. We will see that this solution is unstable. The 
true vacuum has non-zero T and breaks chiral symmetry. 

The potential V^eff (|3.8|) contains a classical term, which is positive and favors T = T = 
0, and a negative one loop term, which becomes larger (i.e. more negative) as \k\ decreases. 
Thus, it is natural to expect that any non-trivial solution of ( |3.9|) will be dominated by 
the low momentum modes of T. 

To find such a solution it is useful to discuss separately two regimes. The first is the 
linear regime, in which 

T{k)T{k) < k'^ , (3.10) 

so that we can expand the log in ( |3.8| ) and keep only the leading quadratic term and hence 
the linear term in the gap equation (|3.9|). In this regime (|3.y|) becomes 



V^ 



(£l±i!)!^r(.) 



+ T(x) = 0, (3.11) 



where V^ is the four-dimensional Euchdean Laplacian. 

As mentioned above, due to the relation with the chiral condensate, ( |3.7D , we expect 
T to depend only on r = vx^. In terms of 

F{r) = (r^ + L2)iT(r) = ^(gl(x) ■ ?i,(0)) , (3.12) 

equation ( p.ll|) takes the form 

F"ir) + -/'{r) + ^;^^^, = ■ (3.13) 

We now argue that the solution of (|3.13|) which is relevant in the linear regime (|3.10|) is 
simply F(r) = C with C a constant to be determined. To see this, note that ( p.l3|) can 
be solved in closed form for r <^ L in terms of the dimensionless coordinate p = r^/kJlJ' 
and for r ^ L in terms of the dimensionless coordinate a = r/A. The solution for r -^L 
which is regular as p ^ is expressed in terms of Bessel functions as 

F{p) - ^-^ , (3.14) 

P 

and approaches a constant for small p. The regime of validity of ( |3.14| ), r <^ L, corresponds 

to p <^ a/A/L ^ 1. Thus, to leading order in A/L, the only part of the solution ( p.l4|) we 

are sensitive to is its value at p = 0. 

For r ^ L the leading behavior of the solution at large a can also be expressed in 

terms of Bessel functions as 

F(a) ~ ^i(HM) , (3.15) 

a 

and approaches a constant for large a. Since r ^ L implies cr ^ 1, we are again only 

sensitive to this constant. 

Thus, F{r) is constant both for r <^ L and for r ^ L. One can show that it does not 
exhibit any non-trivial variation for r r^ L, i.e. the two constants are the same (to leading 
order in A/L). This is a consequence of the fact that the coefficient of F{r) in (|3.13| ) is 
bounded from above by A/L"^, which means that the scale of variation of the solution is 
a/L"^/A, a distance scale that is much larger than L in the weak coupling regime. It can 
also be verified by a numerical solution of ( |3.13| ). 

In terms of T{r) ( |3.12| ), the solution is 

C 
10 



with Fourier transform 

A 2fi -kL 

T{k) = ^^^^ . (3.17) 



Since equation (|3.11|) is linear, the constant C is not determined by the hnear analysis. 
We will determine it below by matching to the non-linear regime. 

As we anticipated from the form of the potential ( p.8|) , the vacuum expectation value 
of T ( p.lTp grows with decreasing momentum. The condition ( |3.1CI| ), which is necessary 
for the validity of the preceeding discussion, is bound to be violated for sufficiently low 
momenta. We will see that this happens in the regime kL <^ 1. Thus, in the transition 
region we can neglect the exponential in (|3.17| ), and find that the linearity assumption 



breaks down at /c ~ y/C. 

We saw before that in the linear regime ( |3.1C1| ) the gap equation ( |3.9|) reduces to ( |3.11|) . 
In the opposite, strongly non-linear regime, 

T{k)T{k) > P , (3.18) 

one instead finds 



.2 I r2 



3 



,4^(x^ + L-j2 .^.^ 



T{k) J A 

The solution of ( |3T9| ) is 

T(x) =T(x) = AS^{x 

T{k) =T{k) = A , 
where the constant A is determined by (|3.19| ), 



rf*x^ — '—T{x)e-"'-'' . (3.19) 



(3.20) 



^ = Vf • (3.21) 



The 5-function behavior of (|3.20|) means that T(r) goes rapidly to zero for large r. 



We can now determine the constant C by matching the linear and non-linear regimes. 
Comparing ( |3.17|) to the second line of ( p.20|) we see that the transition between the two 
occurs around the momentum k* which satisfies 

At that point we expect the parameter 

T{k*)T{k*) A^ 



{k*Y ~ (fc*)2 
11 



1 . (3.23) 



This means that 



^==i£)-4^- (^•2*) 



We are thus led to the picture shown in figure 2. For momenta k ^ yA/L^ or, equivalently, 
distances \x\ <^ a/L^/A, we are in the hnear regime, in which the gap equation takes the 
form (|3.11| ), and its solution is given by ( p.l6| ), ( |3.17| ). On the other hand, for momenta 
k <^ a/A/L^ the system is in the non-linear regime, where the gap equation is given by 
( p.l9|) and its solution by (|3.2C1| ). The momentum scale 

at which the transition between the linear and non-linear regimes takes place is very low: 
k* <^ 1/L for weak coupling, A ^ L. This provides an aposteriori justification for setting 
exp(— /c*L) ~ 1 in the analysis above. 

T(k) 




Fig. 2: Behavior of T(k) as a function of momentum. 



To show that the non-local NJL model (|3.5| ) breaks chiral symmetry, it is important 



to establish that the non-trivial solution of the gap equation constructed above has lower 
energy than the trivial solution T = 0. This can be seen to be a property of any sufficiently 
well behaved solution of the gap equation. 

Indeed, multiplying ( 3.9|) by e^^"^T{y) and integrating over k and y leads to 

i/.W + .^)!TWTW./|l^™)^. (3.26) 

12 



Plugging this into (|3.8|) we find that Ves can be written as 



f d^k \ T{k)T{k) ( T{k)T{k) \ 



(3.27) 



The integrand in ( |3.27| ) is non-positiveQ for any finite T{k). Assuming that the integral 
over k converges, which can be verified to be the case using the explicit form of T(/c) that 
we found, ( p.l7|) , (|3.20|) , one concludes that Ves < for the non-trivial solution of the gap 



equation, i.e. the latter has lower energy than the trivial solution T = 0. 
We finish this section with a few comments: 

(1) It is interesting to use the results above to calculate the quark anti-quark condensate 
{q^ix) ■ qR{y)). Using (|3.12| ), ( |3.16| ) we see that in the linear regime 

{q[{x) ■ qniO)) c^ ^ . (3.28) 

Thus, for \x\ = r '^ 1/k* the chiral condensate is independent of r and A. For 
r ^ 1/k* it goes rapidly to zero. The full solution can in principle be obtained 
numerically to arbitrarily high precision. 

(2) We emphasize again that the results obtained in this section are valid to leading order 
in X/L. One can compute higher order corrections to these results by including the 
non-linear terms in the gauge field Lagrangian and other corrections discussed above. 

(3) The chiral condensate (|3:28|) transforms in the (Nf,Nf) of the U{Nf)L x U{Nf)R 
global symmetry ( p.3|) . The flavor indices are suppressed above. The expectation 
value ( p:28D breaks the chiral symmetry U{Nf)L x U{Nf)R -^ f/(A^/)diag- 

(4) The condensate {q\-qR) that we found in the non-local NJL model exhibits the type of 
behavior one might expect in QCD. There, the chiral condensate should presumably 
be roughly constant for r smaller than 1/Aqcd and go to zero for large r. 

4. Strong coupling analysis 

In the last section we analyzed the dynamics of the chiral fermions qL and q^ in the 
weakly coupled regime A <^ L. As we increase QgNc, or bring the DS-branes closer, the 
effective coupling of the fermions due to exchange of il>4-brane modes increases, and the 
approximations which we employed break down. 



Indeed, defining x = T{k)T{k)/k^, the term in square brackets in ( 3.27 ) is given by [• • •] 
jj^ — log(l + x), which is negative for all x > 0. 
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For L ^ A, there is another weakly coupled description of the dynamics. Instead 
of studying fermions coupled via five-dimensional gluon exchange we should consider the 
dynamics of L>8-branes in the near-horizon geometry of the Nc -D4-branes. In this section 
we will use this description to analyze the question of chiral symmetry breaking in this 
regime. 

Using the conventions of |jT^ in which the radial coordinate U has dimensions of 
energy, the metric is given by 



3 3 



ds' = (^) ' {r^^.dx'^dx'' - (dx'f) - (^] ' {{a'dUf + (a'Ufdnl) , (4.1) 

where O4 labels the angular directions in (56789). The fourbrane geometry also has a 
non-trivial dilaton background, 

The parameter R is given by 

r3 = TxgsN,{a'f'^ = ^N,a' = nXa' . (4.3) 

47r 

In this section we take x'^ to be non-compact. In the next section we will discuss the 
compact case studied in fl^. In what follows we set a' = 1. 



We now consider a probe L>8-brane propagating in the geometry ( |4.1j ).H The D8- 
brane wraps IR ' x 5"^ and forms a curve U = U{x'^) in the (t/, x'^) plane, whose shape 
is determined by solving the equations of motion that follow from the DBI action on the 



-D8-brane. In the background ( [4.1|) , ( f4.2|) the action is 



Sd8 = -TsVs+iV4 I dx''U''\ll+ (^^ [/'2 , (4.4) 

where U' = dU /dx^, V3+1 is the volume of 3 -(- 1 dimensional Minkowski spacetime, and 
V4 the volume of a unit 5"^. Since the integrand of ([4. 41) has no explicit x^ dependence. 



U{x ) satisfies the first order equation [^ 



r/4 

U'o . (4.5) 



i + i^fu- 



The discussion of Nf coincident DS-branes is essentially identical. 
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An interesting solution of (|4.5| ) is a U-shaped curve in the (C/, x^) plane, which approaches 

x^ = ±-j as t/ ^ oo, and is equal to Uq at x^ = 0, the point of closest approach of the 

curve to U = 0. The curve is symmetric under x^ 

Solving (^]^) for U' and integrating leads to 

-u 



4 
-X . 



x\U) 



dU 



^^{if'(^-l 



1/2 



(4.6) 



which can be written in terms of complete and incomplete Beta functions as 



X 



(^) = i^ [^(9/16, 1/2) - B{U^/U'; 9/16, 1/2)] 



(4.7) 



From this we read off the asymptotic value x^(oo) = L/2, 

1-1 1 9 1 
4 ° H6'2^ 



(4. 



The incomplete Beta function has an expansion at small z given by 

B{z]a,b) = z" 



1 1-6 

-H -z + 



a a + 1 
Keeping the first term in this expansion in ( [1.7| ) gives the form of the curve at large U: 



(4.9) 



. 2 R2Ui 



k — rr4 
2 



9^-x 



(4.10) 



The part of the D-brane that corresponds to x^ < is determined by the symmetry 
U{x'^) = U{—x'^). The full Z58-brane is shown in fig. 3. 



(5) 

A 



(6) 



-* (4) 




Fig. 3: A slice of the full DS-hvane configuration showing the D8 and DS-hvane 
joined into a single D8-brane by a wornihole. 
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The solution described above is not the near-horizon description of the brane config- 
uration of figure 1. That configuration corresponds in the geometry (^4.1J ) to a stack of 
Nf L>8-branes stretched in U and locahzed at x^ = —L/2, and a stack of Nf L>8-branes 
stretched in U and locahzed at x^ = L/2. In particular, unlike the solution ( [4.6|) , the D8 
and DS-branes are not connected in this case. Thus, the configuration of figure 1 pre- 
serves a U{Nf)L X U{Nf)fi symmetry acting separately on the D8 and L>8-branes, unlike 
the configuration (j^^ where the two are connected and the symmetry is only a single 
(diagonal) U{Nf). 

It is a dynamical question whether for fixed L it is the separated parallel brane con- 
figuration or the connected curved one, with Uq a function of L ( [4.8|) , that minimizes 
the energy density. The energy per unit volume is infinite, but the relative energy den- 
sity of the two configurations is finite and can be computed by integrating the energy 
difference for each slice dU. After substituting the curved configuration's U' from ( ^.5| ) 
and rewriting the DBI action ( [4.4| ) as an integral over U rather than x^, the difference 
AE = -^straight - EcuTved IS proportioual to 

AE ~ r° (t/5/2 - o) rft/ + r [ t/5/2 - f/5/2 fl - ^^ 

= - -Uy^B{~7/16, 1/2) ^ 0.052 C/q/^ . 

Thus, -^straight > -^'curved, SO the curved configuration is preferred. Physically, this is due 
to the attractive force between the D8 and DS-branes mediated by the L)4-brane fields. 
In the previous section we studied its consequences in the weakly coupled regime. In the 
strongly coupled regime under consideration here, the attractive force leads to a large 
deformation of the DS-hranes which can be seen in figure 3. 

To check the validity of the supergravity approximation, it is convenient to rewrite 
( [4.8| ) in terms of the variables A ( |2.4D and L. Omitting constants of order one, we have 

Uo^^. (4.12) 




As explained in |T^, a necessary condition for the supergravity solution p.l| ) to be valid 
is that the curvature in string units is small, which is the case when the effective 't Hooft 
coupling is large, 

a(^)'»1, (4.13) 
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or, equivalently, 

At/ > 1 . (4.14) 

Uq ( [4 .121) satisfies ( [4.14|) if A ^ L. Fixing L and increasing A puslies Uq furtlier into tlie 
regime of validity of supergravity.El On the other hand, decreasing A leads to smaller Uq, 
and eventually, when A ~ L the curvature at Uq becomes of order one and the supergravity 
description breaks down. 

In section 3 we saw that the weakly coupled non-local NJL description of the four- 
dimensional dynamics is valid for A ^ L. Now we see that the description in terms of 
a curved DS-brane is valid for X ^ L. This is an example of a bulk-boundary duality 
analogous to that between A^ = 4 SYM and supergravity in AdS^ x S^ . The analog of the 
't Hooft coupling of A^ = 4 SYM in our case is the dimensionless ratio X/L. 



As emphasized in [^ , the fact that what looks asymptotically like two disconnected 
stacks of D8 and L>8-branes is in fact part of a connected stack of curved D8-branes 
provides a nice geometric realization of chiral symmetry breaking. At high energies (which 
correspond to large U [|l^) one sees an approximate U^Nf)^ x U{Nf)ji symmetry while 
at low energies (small U) only t/(A^/)diag is manifest. 

The energy scale associated with chiral symmetry breaking in the supergravity regime 
is Uq (|4.12|) . This scale can be thought of as the constituent quark mass of ql, qr. 



Indeed, in the limit X ^ L the spectrum contains free quark states which correspond to 
fundamental strings stretched between the curved DS-branes and U = (the location of 
the Nc il>4-branes). The energy of such strings is Uq. The analog of this scale at weak 
coupling is k* ( |3.25| ). As in other bulk-boundary dualities, the mass goes like a different 
power of the coupling X/L in the weak and strong coupling regimes. 

To make the description of the 3 -|- 1-dimensional dynamics in terms of a curved D8- 
brane in the il>4-brane geometry more precise, we need to find the map between bulk fields 
in the geometry (|4.1|) and boundary operators. Once this map is established, one can use 



the standard tools of holographic dualities to study the boundary theory. In particular, giv- 
ing an expectation value to a non-normalizable operator in the bulk corresponds to adding 
the dual boundary operator to the Lagrangian. Normalizable bulk v.e.v.'s correspond to 
giving expectation values to the dual boundary operators. 



There is an upper bound on Uq coming from the requirement that gs{Uo) «C 1, but it involves 
Nc and we will not discuss it here. 
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In our system, there are two kinds of bulk fields. One is closed string fields (such as 
the dilaton, graviton, etc) in the geometry ( [4.1|) , which exist even when the L>8-branes are 
absent. The other is open string fields on the DS-hranes. Both live in the bulk {i.e. at any 
U) and couple to ql, Qr- We will briefly comment on the bulk-boundary map for some 
open string fields, leaving a more detailed analysis to the future. 

In order to find the boundary operators corresponding to different open string modes 
we go back to the D-brane configuration we started with (figure 1). Consider first the open 
string tachyon stretched between the D8 and DS-branes. This complex scalar field which 
transforms as (A^/, Nf) under U{Nf)LX U{Nf)ji couples to the fermions via a Yukawa-type 
interaction 

CiC:^Tq],-qL + Tql-qR. (4.15) 

Thus, in the near-horizon geometry of the L>4-branes, the open string tachyon T, T is dual 
to the boundary operators 

T ^ q]i-qL , 

_ f (4.16) 

T ^ql-qR . 

In particular, turning on a mass for q^, q^ corresponds to giving a non-normalizable 
expectation value to the field T, while a v.e.v. for gjj ■ q^ corresponds to a normalizable 
expectation value of T. 

As we discussed above, the curved DS-hrane ( [4.6| ) describes a vacuum with non-zero 
expectation value {q]^ ■ Ql)- Therefore, in addition to its curved shape, the D-brane ( |4^ ) 
must have a non-zero normalizable condensate of the 8 — 8 tachyon T. At first sight this 
may seem odd, but in fact something very similar is known to happen for the closely 
related hairpin D-brane [|T7| - p!9| . In our case, the fact that T has an expectation value can 
be seen from ( [4.15| ). Since (gjj ■ g^) is non-zero, there is a tadpole for T localized at the 



intersection. Therefore, it has a non-zero expectation value, which is also localized in the 
vicinity of U = 0. 

Another mode that we can consider is the scalar (^^) ^ that parametrizes the loca- 
tion of the D8-branes in the x^ direction, and its D8-hrane counterpart (X^) . Unlike 
the tachyon, (^X^) j^ ^ do not have Yukawa-type couplings to ql, q^, as can be easily de- 
duced from symmetry considerations. The lowest dimension coupling consistent with the 
symmetries has the form 

C2 ^ [X^)^ {d,Jl) + {X^)^{d,J^) , (4.17) 
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where Jj^, J^ are the U{Nf)L x U{Nf)fi currents. In the free theory at infinite L they are 
conserved, so the coupling (|4.17|) can be neglected. For finite L, chiral symmetry breaking 
imphes that d^J^^ ^ 7^ and proportional to q\ ■ qRq\i ■ Ql (integrated over some of the 
positions). Thus, (|4.17| ) can be written as 

£, - {X'), {q[ ■ im'n ■ 5^)(^,_^) + Mh (^1 ■ ^^^1. ■ ^^)(,,^p (4-18) 

where the subscripts indicate that in the first term q]^ ■ qRq]^ ■ qL are coupled to an adjoint 
of U{Nf)L and a singlet of U{Nf)fi, and similarly for the second term. 

In the near- horizon geometry ( [4.1| ) , ( [4.18|) implies that the bulk-boundary correspon- 
dence is 

. ' (4.19) 

The fact that the operators on the r.h.s. of ( [4.19 ) have non-zero expectation values in 



our brane configuration implies that the scalars [X^^ , (^^) „ must have a non-zero 
normalizable expectation value. This is nothing but the curved shape of the branes, which 
is given asymptotically by ( [4.10|) . It is normalizable since a;^ approaches L/2 as C/ — > 00 



like L/2 - x^ r^U-i. 

It would be interesting to study the correspondence between bulk fields and boundary 
operators in more detail, but we will not pursue that here. We finish this section with 
some comments on the 8 — 8 string, which seems to play an important role in this system. 

In fiat space the ground state of the 8 — 8 string is tachyonic for L less than the string 
length and has positive mass squared for L larger than the string length (as is the case 
here) . One might be tempted to think that it can be decoupled from the dynamics because 
of its string scale mass. However we believe that this is incorrect. 

One way to see this is to study the geometry of a fundamental string stretched between 
the D8 and D8-branes in the near-horizon geometry of the D4-branes. Naively, such a 
string is stretched in the x^ direction at fixed U. Let U* be the fixed value of U and 
L*/2 = x^{U*). Then the coordinate length of the string is L*, and including the warp 
factor in (^7T|), its physical length is 
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L{U) = L* ( ^ ) . (4.20) 



However, this is not the lowest energy configuration. As is familiar from the study of 
warped geometries [2^22], the string dips into the small U region due to the warping. To 



find the amount of this dipping, we can proceed as follows. 
The string is described by the Nambu-Goto action, 

SNG^-Jd'^^V-dethab (4.21) 



where 



hab=daX'''dbX''GMN (4.22) 

is the induced metric on the string and we omitted an overall numerical factor. Plugging 
(ffTl) into the Nambu-Goto action and omitting a factor of the length of time, we find that 




Sng = - I _.t^x^W( o ) +U'^ ■ (4-23) 

The shape U{x'^) satisfies the first order equation 



(S)^ _^f/r^^ 



fl^3..,o \ R 



(4.24) 



i+{§ru- 



Here Uq is the minimal value of U that the fundamental string attains, at x = 0, where 
U' vanishes. It is determined by U* (and L*, but the latter is also determined by U* via 
the shape of the DS-brane, which we take as given) . 

Some algebra leads to the following relation between the various parameters: 

3 ^ 

i* I R V f^ ,-4,3 



2c/r' wi ^1 



(ixx'^ix"^ - 1)"2 . (4.25) 



This equation is valid for all U*, but it simplifies for large U* , for which L* approaches 
( ^TBI ) , and one can assume that U* ^ Uq so that the upper limit of the integral ( ^4.25| ) 
can be taken to infinity. This gives the following relation: 



3 

R •' 



Note in particular that the amount by which the fundamental string descends into the 
il>4-brane throat is independent of U* for large U* . 

We can also compare the minimal value of U for the fundamental string and the 



D8-brane. The latter is given by ( [4.8|) ; combining it with (|4.26| ) we conclude that 

Uo ) - 3i?(9/16, 1/2) - ^■■^** ■ '■*-^^> 

Thus, the fundamental string reaches the vicinity of the tip of the D8-brane, but stays 
away from it by a finite amount. This behavior is shown in figure 4. 
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Fig. 4: Configurations of fundamental strings with ends on tlie DS and D8-branes. 

Another interesting quantity to compute is the length Lcurved of the fundamental string 
whose shape is described by ( [4.24| ) . Dividing this length by the length of the straight string 
(and omitting some numerical constants) one finds that the ratio is 



^curved 



^straight 



w 



u, 



:f) 



(4.28) 



Hence, for large U* the curved string is much shorter than the straight one. 

Thus, even a string that starts at a large value of U and might be expected to naively 
decouple from the infrared physics at small U in fact does not. To minimize its length, the 
fundamental string descends into the £)4-brane throat, reaching the vicinity of the place 
where the L>8-brane itself turns around, and proceeds to the other side in x^. 



5. NJL and QCD 

An interesting generalization of the brane configuration we have been studying is 
obtained by compactifying the x^ direction on a circle of radius i?4, with antiperiodic 
boundary conditions for the fermions living on the L>4-branes. For finite R4 the U{Nc) 
gauge field on the L>4-branes becomes dynamical. The antiperiodic boundary conditions 
give a mass to the adjoint fermions and scalars coming from 4 — 4 strings, at tree level 
and one loop, respectively. Below the mass of the adjoint fermions and scalars the dy- 
namical degrees of freedom are a U{Nc) gauge field and Nf fermions in the fundamental 
representation of the gauge group. 
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This brane configuration was studied in W^ as a description of QCD witfi Nf flavors 



generalizing tlie work of [|2^, wliere tlie system witfiout DS-branes (and tlius witliout fun- 
damentals) was considered as a model for pure Yang-Mills theory. Its low-energy dynamics 
depends on two dimensionless parameters: X/L and L/i?4. In sections 2 - 4 we considered 
the case L/R4 = and studied the dependence of the dynamics on X/L. In this section 
we will qualitatively describe the dependence of the dynamics on L/R4, which varies over 
the range [0, tt]. This will help us to relate the discussion of sections 2 - 4 of this paper to 
that of [|T3[, and to QCD. 

We start with the case where x^ is compactified on a large circle of radius R4 ^ L. 
As in the previous sections, the dynamics depends on the open string coupling QsN^- For 



weak coupling, QgNc <^ 1, the hierarchy of scales is (compare to (|2.5|)) 

A < /s < L < 1^4 . (5.1) 

The fact that R4 is finite means that the classical four-dimensional 't Hooft coupling A4 

is finite as well. In the regime ( |5.1| ), the four-dimensional coupling A4 is small. Dimen- 
sional transmutation generates a dynamical scale Aqcd, at which four dimensional gauge 
interactions become large. 

For small A4, Aqcd is much smaller than the dynamically generated mass of the 
quarks, k* ( |3.25|) . Therefore, the dynamics splits into two essentially decoupled parts. 



The chiral symmetry breaking is still described by the non-local NJL model of section 3. 
The interactions of the quarks with four dimensional gauge fields provide a small correction 
to their properties at the scale k*. In particular, q^ and qn behave as free particles at 
distance scales of order 1/k*. The gauge fields introduce a confining potential for the 
quarks at the much larger distance scale 1/Aqcd- Thus, for finite R4 we expect to find 
a discrete spectrum of meson resonances with energies of order 2/c* and mass splittings of 
order Aqcd- 

As the open string coupling QsNc increases, the five-dimensional coupling and Aqcd 
increase as well. For 

Zs < L < A < i?4 (5.3) 



the four-dimensional coupling ( |5.2|) is still small, but the five dimensional one is large. 
Thus, one can use the supergravity description of section 4 (with small corrections due to 
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the finiteness of R4) to describe the dynamical mass generation of the quarks, while the 
confining interactions of the quarks with the gauge field are still described by Yang-Mills 
theory with a small 't Hooft coupling ( |5.2|) . In particular, the dynamics still splits into 
the chiral symmetry breaking part and the confining part, which occur at different energy 
scales. 

Further increasing QsNc we reach the regime 

/« < i^ < i?4 < A . (5.4) 

Here, both the four- dimensional and the five-dimensional 't Hooft couplings are large, and 
one needs to use supergravity to study both chiral symmetry breaking and confinement. 
The near-horizon metric of the L>4-branes in this regime is given by 

ds^ = (^) ' {v^^dx^dx'' - f{U){dx^)^) - 
where 



f{u) = i '■^- 



u 

R is given by ( [4 .31) and 



a'uy 

R ) 


(5.5) 


K ^ 

1 


(5.6) 



^-- = TS! • f^-^' 
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As is familiar from studies of near-extremal £)4-branes, the geometry ( p.5| ) is smooth. The 
U coordinate is restricted to t/ > Ukki such that the warp factor in IR ' as well as the 
volume of the four-spheres remain finite for all V . When (|5.4| ) is valid, both [/q ([4.12| ), and 
Vkk (^]3) ^^^ ^^ ^^ region where the curvature of the metric (|5.5| ) is very small, so one 
can study low-lying open and closed strings using supergravity. 

Equation ( p^) implies that Uq ^ Ukk- Therefore, the dynamics associated with 
chiral symmetry breaking is still insensitive to the presence of dynamical four dimensional 
gauge fields. Indeed, in studying the shape of the L>8-branes in section 4 we took f{U) = 1, 
rather than ( ^.6|) . However, this shape is only sensitive to f{U > Uq)- Thus, the analysis 
of section 4 is valid for finite Ukk, up to small corrections that go like {Ukk/UqY . In 
particular, the dynamically generated mass of the quarks is still given by Uq (^T^), ( [4.12| ). 

What does change significantly in this case is the properties of constituent quarks, 
which for Ukk = were described by fundamental strings going from Uq down to t/ = 0. 
Indeed, consider a quark anti-quark pair, which corresponds to a fundamental string whose 
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ends lie on the DS-branes, and are separated by a distance a; in IR . For i?4 ^ x ^ L, this 
looks like a string going from the DS-brane down towards U = and another one with 
opposite orientation a distance x apart. Its shape is insensitive to the finite value of R4. 
This implies that the quark and anti-quark are weakly interacting at these scales. 

When the distance x becomes of order R^, the string notices that it can descend no 
further in U and its properties change. This corresponds to the confining potential created 
by the four dimensional dynamics at the distance scale i?4. Thus, in this regime the quarks 
can be thought of as particles with the mass Uq ( [4.12|) which are free at short distances 
and form bound states whose size is of order i?4. 

To summarize we find that if we compactify the system discussed in sections 2 - 
4 on a large circle of radius R4 ^ L, the chiral symmetry breaking dynamics of the 
fermions ql and qr remains essentially unchanged and is almost decoupled from that of 
the four dimensional gluons of U{Nc). At the energy scales associated with chiral symmetry 
breaking the four-dimensional coupling is very small and one can neglect the dynamics of 
the gauge fields. The scale Aqcd at which the gauge coupling becomes strong is well below 
the dynamically generated masses of the quarks. Gauge dynamics leads to the formation 
of bound states whose typical size is 1/Aqcd- 

To get a theory that looks more like QCD we need to take the parameter L/R4 to 



be of order one. For example, [0 discussed the case in which the D8 and L>8-branes are 
maximally separated on the circle, which corresponds to L = nR^. When L and R4 are 
comparable, the only dimensionless parameter that we can vary is X/L ~ A/-R4. 

For X <^ L, the low-energy theory on the branes is QCD with Nf massless fundamen- 
tals, and a small four- dimensional 't Hooft coupling (|5.2|) . The QCD scale is well below 
the other scales, 1/L, I/R4. This is the theory one would like to solve, but unfortunately, 
there are no good analytical tools to study it in this regime. 

For X ^ L the theory is not quite QCD, but it is expected to be in the same universal- 
ity class. In particular, qualitative phenomena such as dynamical symmetry breaking, and 
the existence of a discrete spectrum of massive glueballs and mesons should not change as 
one varies X/L, although the details of the spectrum may change. The limit X ^ L can 
be studied using the supergravity description ( |5.5D . This was done in [|T^,^ and we will 
not review the details here. 

The only point we would like to mention is that unlike the limit L/R4 -^ discussed in 
section 4, in the smooth spacetime ( p.5| ) a configuration with separate D8 and L>8-branes 
terminating at the origin does not make sense, and only the "hairpin" shape joining the 
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asymptotic D8 and Z)8-branes is possible. Thus, in contrast to the discussion of section 
4, where chiral symmetry breaking was simply favored energetically, here it is mandatory. 
This is perhaps not surprising given that the 't Hooft matching conditions require chiral 
symmetry breaking for specific values of Nf [^5|] . 

The discussion of this section leads us to one of the main points of this paper. In 
the case without quarks studied in ||5^, the dynamics depends only on one dimensionless 
parameter, A/i?4. For small values of this parameter one finds pure Yang-Mills QCD, 
which is hard to treat analytically, while for large values one finds a system that differs 
from QCD but can be analyzed using supergravity 

In the -D4 — D8 — D8 system, which contains massless quarks, there are two dimen- 



sionless parameters, X/L and L/R4. For L/R4 ~ 1, the case discussed in [ITsiM , the 
situation is as in ||2^, but for small L/R4 one gets a theory which is solvable both at 
weak coupling X <^ L and at strong coupling A ^ L, by using the non-local NJL model of 
section 3 and the supergravity analysis of section 4, respectively. Moreover, we presented 
evidence that the system is in the same universality class (or phase) for all values of these 
parameters. Thus, we see that the brane construction interpolates between a regime in 
which the dynamics of the quarks is described by the non-local NJL model, and one where 
it is described by QCD. This might help explain why the NJL model is useful in studying 
mesons in QCD [|-|. 

6. Discussion 

Much of the early work describing large Nc QCD using string theory has focused on 
pure Yang-Mills theory, without quarks. In such theories one can compute the spectrum 
of glueball states, but the results must then be compared to lattice QCD since glueballs 
have not been unambiguously identified experimentally. In trying to obtain more realistic 
models of QCD it is important to incorporate quarks and their meson bound states. The 
D4 — D8 — D8 model studied in [p!3| , |2^ and in this paper is a natural construction which 
incorporates the chiral symmetry of QCD with massless flavors, describes the spontaneous 
breaking of this symmetry, and leads to a spectrum of meson bound states. 

Previous work on this model fl^,^ focused on the region in the parameter space of 
brane configurations in which the separation between the L>8-branes, L, is of the order of 
the radius of the extra dimension along the Z)4-branes, R4. If the five-dimensional 't Hooft 
coupling A is small, A ^ L.R^, this brane configuration describes QCD, and is difficult to 
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analyze. For large A the dynamics is not quite that of QCD but it can be analyzed using 
supergravity. There are reasons to believe that the theories one gets in the two limits are 
in the same universality class. 



In this paper we studied the brane configuration of [1T3| , |2^ in the opposite limit R4 ^ 
L. Consideration of this limit and the space of brane configurations in general leads to 
some new insights into the dynamics of quarks and mesons. 

One of the surprising results of our analysis is that there exist four-dimensional models 
of quarks without dynamical gauge fields which exhibit non-trivial infrared dynamics. In 
local quantum field theory this is believed to be impossible. Our model contains a non- 
local interaction between the left and right-handed quarks, but since this interaction arises 
in a D-brane system, we are assured that it does not lead to any pathologies. 

The model in question is a non-local analog of the Nambu-Jona-Lasinio model, which 
exhibits dynamical symmetry breaking for arbitrarily small value of the coupling. While 
we obtained this model from string theory, one can construct it directly in field theory as 
follows. Consider a five- dimensional U{Nc) gauge theory with two codimension one defects 
separated by a distance L. At the two defects there are Nf left and right-handed fermions 
in the fundamental representation of the gauge group, ql and qr. The dynamics of the 
four-dimensional fermions and five-dimensional gauge field is governed by the action ( |3.4|) . 
The model has a U{Nf)L x U{Nf)ii global symmetry. Moreover, since the U{Nc) gauge 
fields are higher dimensional, U{Nc) acts as a global symmetry as well. 

The dynamical degrees of freedom in the four-dimensional theory are the fermions 
ql and qr. The five-dimensional gauge field provides an effective coupling between them. 
This coupling becomes weaker as the distance L increases, since the five dimensional gauge 
theory is infrared free. In the analogy to superconductivity that motivated 0, ql and qr 
are analogous to the electrons, and the five-dimensional U{Nc) gauge bosons are analogous 
to the phonons. 

To study the interaction among the fermions it is convenient to integrate out the five- 
dimensional gauge field. In the limit L 3> A, one can do that in the leading, single gluon 
exchange approximation. This gives rise to the non-local NJL model ( |3.5| ), which breaks 
chiral symmetry for arbitrarily large L [i.e. arbitrarily weak coupling). 

Since the non-local NJL model ( |3.5| ) contains only vector degrees of freedom of U{Nc), 
it is exactly solvable in this limit. It can be thought of as a generalization to four dimensions 



of models such as the Gross- Neveu model [Q, and the 't Hooft model [|2^ of two-dimensional 
QCD. 
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Compactifying x"^ on a circle of radius i?4 provides a continuous interpolation between 
the NJL model ( p.5| ) and QCD. The latter is obtained in the limit i?4 ~ L ^ A. Our 
analysis of the NJL model strongly suggests that it is in the same universality class as 
QCD. 

Physically, the reason that the above construction simplifies the study of quark dy- 
namics is the following. The basic mechanism for chiral symmetry breaking and the binding 
of quarks into mesons is the attractive interaction between the quarks due to exchange of 
U{Nc) gauge bosons. Our model allows one to separate the two scales associated with 
chiral symmetry breaking and confinement. In the NJL limit, the mass scale associated 
with chiral symmetry breaking is much higher than that of confinement, and all the com- 
plications associated with the latter disappear in studying the former. In the QCD limit 
the two scales are comparable, which makes the model harder to analyze. 

Many interesting questions remain to be addressed. As in previous discussions, we 
have dealt with the theory in the limit of zero bare quark mass. Since one can clearly 
make sense of the quark mass perturbation in the NJL model, we expect the same to be 
true on the supergravity side. It would be interesting to identify this perturbation and to 
study the theory with massive fiavors in the different regimes. This is likely to involve a 
better understanding of the 8 — 8 string discussed in section 4. The detailed correspondence 
between bulk and boundary fields also remains to be worked out. 

The NJL model has been widely used to study a host of phenomenological issues in 
QCD which are not amenable to perturbation theory, including the behavior of the theory 
at finite temperature and quark density. It would be interesting to address these issues in 
the non-local (gauged) NJL model that arises here as well as in its string theory dual. 
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Appendix A. Conventions and useful results 

We use a "mainly minus" metric convention. In particular, the four-dimensional flat 



space metric is 77^,^ = diag(l, — 1, — 1, — 1). We use the same conventions as [^ for the 
Weyl representation of gamma matrices in four dimensions. Thus 

^-(^. ;) (-) 

with a'^ = (1, a) and a'" = (1, -a). 

We require two Fierz identities which rely on the two relations 

(A. 2) 



the flrst of which is eqn (3.77) of [^^ (which also holds with a'^ replaced by a^) and the 
second follows from the flrst using eqn (3.80) of |27]. 



Using these, we have for Grassman-valued Weyl spinors 

(V'Il^''^2l) (^L^M^4i?) = -2 (i^lLi^^R) (^LV'2l) (A.3) 

and 

Fourier transforms in d spacetime dimensions are given by 



f{k) = I d''xe'^-''f{x) , 
m = I (0ze---/(fc) 



(A.5) 



Five-dimensional Fourier transforms of spherically symmetric functions of k in Euclidean 
space, 

^^i^d^De^'-^ , (A.6) 

can be computed using spherical coordinates to give 

7^/ ^ 1 /"°° „ , 4~., , /sin/cx cos/cx\ ,. _, 

f (,) = _ y^ rftt^f (,) [-— - __ j . (A.7) 

As an example, if F{}z) = l//c^ we flnd the coordinate space propagator 
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